Decoherence of flux qubits due to 1// flux noise 
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We have investigated decoherence in Josephson-j unction flux qubits. Based on the measurements 
of decoherence at various bias conditions, we discriminate contributions of different noise sources. 
We present a Gaussian decay function extracted from the echo signal as evidence of dephasing 
due to 1/f flux noise whose spectral density is evaluated to be about (10~^$o)'^/Hz at 1 Hz. We 
also demonstrate that at an optimal bias condition where the noise sources are well decoupled the 
coherence observed in the echo measurement is mainly limited by energy relaxation of the qubit. 
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A Josephson-junction qubit is a macroscopic quantum 
object realized in superconducting electric circuits Q. 
As a solid-state device sitting in the midst of the circuit 
and on top of the substrate, the qubit inherently inter- 
acts with various degrees of freedom which decohere the 
qubit. For understanding and suppressing the decoher- 
ence mechanisms, it is important to identify the domi- 
nant noise sources among the many possibilities. 

A strategy to tackle this problem is to use the tun- 
ability of the qubit. By changing a bias parameter, one 
can control the qubit energy as well as the coupling be- 
tween the qubit and a particular degree of freedom in the 
environment Therefore, the measurement of the de- 
coherence as a function of the bias parameters provides 
valuable information about the environmental degrees of 
freedom. Such a strategy has been applied to charge- 
phase qubits 0, Q, charge qubits 0, p , flux qubits Q , 
and phase qubits . 

In this Letter, we report a detailed characterization of 
decoherence in flux qubits based on a similar approach. 
In a range of bias conditions, we identify flux noise as the 
dominant source of the qubit dephasing. We observe de- 
phasing due to 1/f flux noise which gives rise to a charac- 
teristic Gaussian decay function that seems to be hidden 
in the previous studies by other decoherence mechanisms 
0, . We also quantify the noise spectrum density. On 
the other hand, we show that at an optimal bias condi- 
tion, where those noises are effectively decoupled from 
the qubit, the echo decay rate (~ (3.5 ns)~^ in one sam- 
ple) is about a half of the energy relaxation rate, and is 
mainly limited by the latter process. 

The effective Hamiltonian of a flux qubit can be writ- 
ten as Hq = —^Uz — ^'^x^ where cr^, are Pauli matrices 
[lo| . The energy difference between the two eigenstates is 
given by E^i = Ve^ -I- A^. Here A is the tunnel splitting 
between two states with opposite directions of persistent 
current along the loop, e — 2/p$o(n0 — n*^) is the en- 
ergy bias between the two states, and Ip is the persistent 
current along the qubit loop. We define the normalized 



magnetic flux in the qubit loop rt0 = 'I'cx/^o, using the 
externally applied flux through the qubit $cx and the 
flux quantum <f>o. As the flux qubit we study is cou- 
pled to a SQUID for readout (Fig. P), the SQUID bias 
current also affects flux in the qubit. The normalized 
/f,-dependent flux is thus defined as n*^{Ib) = 0.5— $/j,/$o 
so that 71,0 = n*^ corresponds to e = 0, i.e., degeneracy 
of the two persistent current states. Note that affects 
the qubit only through the fiux <I>/j,(/f,) which is induced 
by the circulating current in the SQUID loop 

We characterize the coherence of the qubit through 
measurements of energy relaxation, free-induction decay 
(FID), and echo decay. From those measurements, we 
deduce the energy-relaxation rate as well as the pure de- 
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FIG. 1: Scanning electron micrograph of a sample and a 
sketch of the measurement setup. The qubit shares a part 
of the loop with a SQUID shunted by an on-chip capaci- 
tor. A current bias line (lb) and a voltage measurement line 
(Vm) consist of lossy coaxial cables and are connected to the 
SQUID via on-chip resistors. Microwave current pulses are 
fed through an on-chip control line (Imw) inductively cou- 
pled to the qubit. The sample chip is enclosed in a copper 
shield box, and all the wires are electrically shielded as well. 
CPF stands for a copper-powder filter. Magnetic flux bias is 
applied with an external superconducting coil connected to 
a battery-powered current source. The sample is cooled to 
20 mK in a dilution refrigerator magnetically shielded with 
three ^-metal layers at room temperature. 
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FIG. 2: (color online). Measurements of (a) energy relax- 
ation, (b) free-induction decay, and (c) echo decay at the op- 
timal bias condition — nj, and It — ■ Red plots are 
for sample A, and black ones are for sample B. The insets 
depict the applied pulse sequence for each measurement. In 
the measurement of free-induction decay, the microwave car- 
rier frequency was detuned from the qubit frequency Eoi/h 
by 10 MHz. Solid curves in (a) and (c) are exponential fits. 



phasing components in the decoherence. Pure dephasing 
is caused by the fluctuations of Eqi which are induced 
by the fluctuations of parameters A, in the Hamiltonian. 
Here 's could be ti^ and /(, as well as other implicit pa- 
rameters inside A and /p, such as the island charge offsets 
and the junction critical currents in the qubit. Through 
the two externally controlled parameters and lb, we 
vary the coupling factor and observe the change in 
the dephasing. Since the derivatives depend differ- 
ently from each other on the control parameters and 
/{,, and assuming that the fluctuations of A^ are inde- 
pendent of the control parameters, we can separate the 
contributions in the dephasing from different sources. 

The experiments were done with samples fabricated by 
electron-beam lithography and shadow evaporation of Al 
films (Fig. nj. The qubit is the small superconducting 
loop intersected by four Josephson junctions [l^ . among 
which one is smaller than others by a factor of ^^0.56. 
We have studied five similar samples, and here data from 
two representative samples are presented. Spectroscopic 
measurements revealed that Ip and A/h were 0.16 fiA 
and 5.445 GHz for sample A, and 0.34 ^A and 5.08 GHz 
for sample B, respectively. The SQUIDs had maximum 
critical currents of ^^3.6 /lA and ^7 /iA, and the switching 
currents under the flux bias used in the qubit readout 
were ~1.6 fxA and ~4.6 fxA, respectively. 

The qubit is initialized to the ground state by relax- 
ation, and is controlled by a sequence of resonant mi- 
crowave pulses. The amplitude and the width of the 
pulses are adjusted so that the desired qubit rotations 
are implemented. At the end of the pulse sequence, the 
qubit state is read out by using the SQUID lia. As the 
qubit in its eigenstates has a finite expectation value of 



the persistent current (except for = np, critical cur- 
rent of the SQUID is slightly varied corresponding to each 
eigenstate. Accordingly, the switching probability of the 
SQUID into the voltage state, responding to a current 
pulse with a given amplitude and a width, is different for 
each eigenstate of the qubit Q . The switching probabil- 
ity Psw, which thus reflects the population of the qubit 
eigenstates, is obtained by repeating the procedure — 
initialization, control, and readout — a number of times 
(typically 10*) and counting the switching events. 

Figurc[21shows the results at the optimal bias condition 
71,0 = 71 J and Ih ~ I^, where both g^"'- and were set 
to zero so that the dephasing due to fluctuations of 
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and /ft is minimal 0, Il2 | . Wc find that energy- relaxation 
and echo signals decay exponentially with rates Fi and 
T2E, respectively, and moreover T2E — ri/2. Therefore, 
at this bias condition the coherence observed in the echo 
measurement is mainly limited by the energy-relaxation 
process. In the analysis of Psw as a function of the read- 
out current pulse height (data not shown), wc have also 
found that the qubit mostly relaxes to the ground state, 
not to the mixture of the ground and excited states. This 
excludes any classical noise from the possible source in- 
ducing the relaxation, and suggests that the spontaneous 
decay of the qubit is the dominant process. 

Next, we study how the fluctuations of /{, affect the 
coherence of the qubit. Wc show that we can dis- 
criminate the contribution from that of 7i0 fluctuations, 
even though the /(, fluctuations also couple to the qubit 
through the flux degree of freedom ti^. First, we present 
the result of qubit spectroscopy at various /;,. For each 
III we swept 7i0 and found a minimum of Eqi to deter- 
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FIG. 3; (color online). Decoherence at various SQUID bias 
currents It (sample A). Qubit flux bias is kept at = n*^. 
(a) Spectroscopically determined vs Solid curve is a 
polynomial fit. (b) Energy-relaxation rate Fi vs (c) De- 
phasing component (e*'^'*')_B in the echo measurements for 
different values of lb- Solid curve is an exponential decay, 
exp(— r^et). (d) FID dephasing rate F<pF, and echo dephas- 
ing rate F^b as a function of /;,. See text for the explanations 
of the fitting curves in (b) and (d). 
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mine which is plotted in Fig. Ola). Notice that at 
= n^, Eoi = A and = 0. The plot also rep- 

resents how Ih affects the ffux in the qubit and how the 
fluctuations of arc transformed into the fluctuations of 
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via -HT^. The interpolated line n*Alb) has a minimum 



at = 0.156 /^A. At If, = the noise current coming 
through the SQUID bias line is split symmetrically into 
two branches of the SQUID and does not produce any 
fluctuations of the flux in the qubit, i.e., = 0. The 
small offset of from zero can be attributed to a small 
asymmetry in the SQUID junction parameters. 

We find a few characteristic behaviors in the decay 
measurements at the bias conditions along the line ~ 
n*^{Ib) (Fig- El- (i) The energy-relaxation and echo de- 
cay curves are always exponential for different values of 
lb, while the rates Fi and T2E largely depend on (h) 
The ratio of pure dephasing rates in the echo measure- 
ment {T^E = F2B — Fi/2) and the FID measurement 
{TipF = F2F — ri/2; with F2_f obtained from exponen- 
tial fitting of FID envelope ^3)' ^fE/^ipF, is signifi- 
cantly smaller than 1 at /f, w while it approaches 
1 at larger |/b — |. These two facts, the exponential 
dephasing and the inefficient recovery of the coherence 
in echo measurement, suggest that the fluctuations in- 
troduced by lb have a relatively uniform spectrum at 
least up to r2£;/(27r) --10 MHz. (iu) The /^-dependent 
part of Fi and F^b are nicely fltted with a term pro- 
portional to (^7^)^ and (^7^)^, respectively. This im- 
plies that the J;,-dependent contributions are solely due 
to lb fluctuations and that the fluctuations are well de- 
coupled at I^: The Fi can be compared with the expec- 
tation from Fermi's golden rule at ~ nj, i.e., Fi = 

M-§^?Cii?Si,[^) M Therefore, we can esti- 

mate the current noise spectrum at the qubit frequency 
Si^{^ = 27r X 5.445 GHz) = 4.4 x 10"^^ A^s/rad corre- 
sponding to the zero-point fluctuations of an impedance 
^ 270 f2. Also, the fourth power term in T^e can be un- 
derstood as a result of the two-fold quadratic couplings 
between Eqi and n^, and and h [l3.ll7| . 

Having seen that the lb fluctuations can be well decou- 
pled at lb =11, we now focus on the effect of fluctu- 
ations. Compared to the measurements above, measure- 
ments at ^ n*^ show qualitatively different behaviors 
(Fig. 01): (i) The energy-relaxation rate Fi is almost con- 
stant in the range of n^, while the dephasing is strongly 
enhanced as \ntj, — n*^\ is increased a little, (ii) The 
echo decay curves are not purely exponential any more, 
and the extracted dephasing component (e"''*^*^)E is flt- 
ted with a Gaussian decay, exp[— (F^^^ t)^] (iii) The 
decay rate F^^ increases almost linearly with — ^^^l- 

Because of the strong n^-dependence of F^^ at 
around n^, we can rule out possible contributions of 
charge and critical-current fluctuations in the qubit. 
Indeed, the observed behaviors are fully consistent 
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FIG. 4; (color online). Decoherence at various flux biases n^i, . 
SQUID bias current is kept at It — I^. (a) Dephasing compo- 
nent {e'*'''')_B in the echo measurements for different values of 
710 (sample A). Solid curve is a Gaussian fit, exp[— (F^^ f)^]. 
(b) Energy relaxation rate Fi and echo dephasing rate F 
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Solid line is a | | fit to F^^^. (c) FID dephasing rate 
Solid line is a \^^\ fit. (d)-(f) Data for sample B. 



with the presence of 1// flux noise. Let us as- 
sume Gaussian fluctuations of with the spec- 

trum density S'„^(a;) = TIjf' Then, we calculate 

/piv{t)\^ _ „™r (dEoi\2 foo J, , c / A sin^(^t/4) i _ 



exp[-(F^^i)'], where F^^ = 



In 2 



dEpi 



0. The 



Gaussian decay function is characteristic of the 1// noise 
spectrum which has a singularity at low frequency. The 
n^-dependence of F^^ follows | | = , which is 

basically linear in the range of studied, as shown in 
Figs. 4(b) and (e). 

Free-induction decay measurement also supports the 
interpretation. As shown in Figs.^Jc) and (f), the Gaus- 
sian decay rate F^^^ also increases almost linearly with 
|7i0 — njl The asymptotic value of the ratio F^^/F^^ 
at large \n^ — n'^\ was between 4.5 and 7.5 for all the sam- 
ples measured, which demonstrates that the echo proce- 
dure largely suppresses the dephasing and that the flux 
noise has a signiflcant amount of low-frequency spectral 
weight. The ratio F^^/F^^ was slightly larger than ex- 
pected {^\n^J\n2 - 5 H) for 1/f noise with an in- 
frared cutoff frequency ujir/{2T:) ^ 1 Hz determined by 
the measurement procedure. We could imagine there ex- 
ist some additional low-frequency fluctuations contribut- 
ing mainly to F^^^. An alternative scenario is an ultravio- 
let cutoff which would make the ratio larger, as suggested 
in Ref.3 for 1/f charge noise. However, such a cutoff 
should modify the Gaussian decay function of (e"^'^*')^;, 
which we did not observe. 
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We evaluated the amplitude of the 1// flux noise from 
the lH^I fitting of T^^, and obtained = 9.3x 10"'^ 

and 9.9 x 10~^ for samples A and B, respectively [l^ . 
Three other samples also gave (1.1-2.0) x 10~^ repro- 
ducively, though the origin of the 1// flux noise remains 
to be understood. One may suspect that the critical- 
current fluctuations of the SQUID junctions produce 
the 1// flux noise via the change of the SQUID circu- 
lating current. However, the normalized 1// critical- 
current fluctuations of the reported value ~ (10^®)^/Hz 
at 1 Hz 113, 113 would give at most 1% of the ^A„^ 
for the mutual inductance of ^20 pH between the qubit 
and the SQUID. It may be worth mentioning that in 
SQUIDs with much larger dimensions similar amount of 
low-frequency flux noise has been known for a long time 
[2^ and not yet understood either. 

It is also interesting to compare the present result with 
a previous work on a flux qubit where the dominant de- 
phasing was attributed to thermal fluctuations of pho- 
ton numbers in the plasma mode of the readout SQUID 
0. There are a few noticeable differences in the result 
of Ref.0: (i) The echo decay was exponential, (ii) r^^^; 
scaled with (^^)^ as a function of n^. (ni) The echo 
decay rate was smallest at 7^ when 7^ In 
our experiment. was smallest at = for any 
/{, (data not shown). Indeed, our estimation shows that 
even in the presence of the same amount of 1 // flux noise 
as in our samples, the contribution would have been con- 
cealed by the photon-number fluctuations in the sam- 
ple in Ref.j3|. In the present work, on the other hand, 
the dephasing rate due to the photon-number fluctua- 
tions is more than an order of magnitude smaller. This 

is mainly because the coupling factor which in- 

creases rapidly when the bias flux through the readout 
SQUID approaches half-integer flux quanta, is about five 
times smaller in our samples. 

In conclusion, we have studied decohcrencc in flux 
qubits. We presented qualitative difference between de- 
cay functions at different bias conditions and charac- 
terized the noises contributing to the decoherence. We 
demonstrated that 1// flux noise constituted one of the 
main decoherence sources. On the other hand, the effect 
was suppressed at the optimal flux bias point where the 
energy relaxation process turned out to be the limiting 
factor of the coherence in the echo decay. Reducing the 
energy relaxation rate is the key issue in improving the 
coherence of flux qubits further. 
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